Prospects and Blind Spots for 
Neutralino Dark Matter 



415 



(N 

O Clifford Cheung^ Lawrence J. Hall^, David Pinner^, and Joshua T. Ruderman^ 
(N 
> 

O ^California Institute of Technology, Pasadena, CA 91125 

^ '^Berkeley Center for Theoretical Physics, University of California, and 

^ Theoretical Physics Group, Lawrence Berkeley National Lab, Berkeley, CA 94720 

<N 

Abstract 

Oh 

P_| Using a simplified model framework, we assess observational limits and discovery prospects 

O for neutralino dark matter, taken here to be a general admixture of bino, wino, and Higgsino. 
Experimental constraints can be weakened or even nullified in regions of parameter space near 
^ 1) purity limits, where the dark matter is mostly bino, wino, or Higgsino, or 2) blind spots, 
^ where the relevant couplings of dark matter to the Z or Higgs bosons vanish identically. We 
^ analytically identify all blind spots relevant to spin-independent and spin-dcpcndcnt scattering 
1^ and show that they arise for diverse choices of relative signs among Mi, M2, and /x. At present, 
XENONIOO and IceCube still permit large swaths of viable parameter space, including the 
^ well-tempered neutralino. On the other hand, upcoming experiments should have sufficient 
^ reach to discover dark matter in much of the remaining parameter space. Our results are 
^ — I broadly applicable, and account for a variety of thermal and non-thermal cosmological histories, 
> including scenarios in which neutralinos are just a component of the observed dark matter 
today. Because this analysis is indifferent to the fine-tuning of electroweak symmetry breaking, 
our findings also hold for many models of neutralino dark matter in the MSSM, NMSSM, and 
Split Supersymmetry. We have identified parameter regions at low tan /3 which sit in a double 
blind spot for both spin-independent and spin-dependent scattering. Interestingly, these low 
tan/3 regions are independently favored in the NMSSM and models of Split Supersymmetry 
which accommodate a Higgs mass near 125 GeV. 
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1 Introduction 



In many supersymmetric theories, the hghtest superpartner (LSP) is a stable neutrahno which 
can account for the observed dark matter (DM) in the Universe. However, two complementary 
experimental efforts seemingly cast doubt on this possibility, at least for the simplest case where 
the LSP is a linear combination of the bino, wino, and Higgsino, x ~ {^iW., h). 

First, a variety of searches from the Large Hadron Collider (LHC) have placed compelling 
limits on supersymmetry, with constraints especially stringent when the LSP is a stable, weakly 
interacting massive particle (WIMP). Naively, this casts doubt on the so-called "WIMP miracle." 
However, the powerful null results from the LHC apply to colored superpartners decaying to 
missing energy — direct limits on the bino, wino, and Higgsino still remain weak. So while the 
LHC challenges supersymmetry as a solution to the hierarchy problem, it does not impose 
strong, direct, constraints on the origin of DM. 

Secondly, after decades of improving technologies and increasing target masses, the direct 
detection of galactic DM has reached unprecedented levels of sensitivity, as shown in Fig. (1). At 
present, the best limits on the scattering of DM against target nuclei are from the XENONIOO 
experiment [1], which probes spin- independent (SI) and spin-dependent (SD) scattering. Com- 
plementary and in some cases more powerful constraints on SD scattering have also been ob- 
tained by the IceCube observatory [2] , which searches for high energy neutrino signals originating 
from DM accumulating inside the sun. In many theories — for example neutralino DM — SI and 
SD scattering is mediated by Higgs and Z boson exchange, respectively. Cross-sections corre- 
sponding to different values of the couplings Ch^x ^^id Czxx shown in Fig. (1), which reflect 
the fact that asi oc c\^^ and cjsd oc c|^^. For neutralino DM, both couplings originate from 
the electroweak gauge couplings g' ~ 0.35 and g ~ 0.65, so the naive conclusion of Fig. (1) is 
that neutralino DM is presently excluded. However, this argument against neutralino DM is 
incorrect. The couplings Chxx ^^"^ ^Zxx ^-^ise from multiple contributions of the same order which 
can constructively or destructively interfere; upon squaring the resultant couplings, one finds 
SI and SD cross-sections which can easily be enhanced or suppressed by an order of magnitude 
from the naive expectation. While current constraints do not require any particular fine-tuning 
among parameters — ^just relative signs — this will not be so easy for future limits. 

In this paper, we explore the observational status and future of neutralino DM in the context 
of simplified models, defined here to be minimal theories of weak scale SUSY, decoupling all 
but a handful of superpartners relevant for DM phenomenology. Crucially, these models are 
characterized by a small number of theory parameters defined at the weak scale. Our aim is to 
identify the regions of parameter space that are presently allowed and understand the detection 
prospects for upcoming experiments. 

A priori, the parameter space of supersymmetric DM is vast, which is why typical analyses 
employ exhaustive parameter scans. However, the majority of these parameters are irrelevant to 
neutralino DM if we assume that the scalar states — heavy Higgs bosons, squarks, and sleptons — 
are sufficiently decoupled to not significantly affect processes relevant to the cosmological or 
observational properties of DM. This simplified model approach is further motivated by the 
absence so far of signals at the LHC for supersymmetry and heavy Higgs states. Similarly, LHC 
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data motivate the assumption that the SI scattering of DM against the nucleon is dominated 
at tree-level by the exchange of a standard model (SM) Higgs of mass near 125 GeV. These 
assumptions lead to a manageable parameter space for neutralino DM which an admixture of 
bino, wino, and Higgsino. This parameter space is comprised of mass parameters for b, w and h 
and the ratio of vacuum expectation values, 

X-^{b,w,h) : (Mi,M2,/i,tan/3). (1) 

We find that our analysis is little altered for squarks as light as 1 TeV, and we quantify their 
effects for lower mass values. For simplicity, throughout this work we assume CP conservation 
and we consider > itt-w^ so that DM is above threshold to annihilate into W^W~ . In 
particular, we do not consider the Z/h pole regions, and we do not consider light DM, < 
10 GeV. 

This minimal framework for neutralino DM is relevant for a remarkably wide range of the- 
ories. The most obvious case is the minimal supersymmetric standard model (MSSM) with 
large tan fi and multi-TeV squarks to yield a 125 GeV Higgs boson. On the other hand, x can 
also describe neutralino DM in the next-to- minimal supersymmetric standard model (NMSSM), 
provided that the singlino component of x is sufficiently small. Such theories can account for a 
125 GeV Higgs boson with less fine-tuning in electroweak symmetry breaking than the MSSM, 
and prefer smaller values of tan /3 [8] . Our setup also characterizes DM in a variety of theories of 
Split Supersymmetry [9] , with tan (5 typically decreasing as the mass of the scalar superpartners 
is raised. In Split Supersymmetry, the decoupling of the scalars is guaranteed, while in other 
theories this is a greater assumption. Notably, the existence of light DM does not invoke any 
additional tuning, since the neutralino mass is protected by chiral symmetries. In this paper we 
are agnostic about the role of supersymmetry for the naturalness puzzle — as such, our analysis 
applies to all of the above theories. 

Our main conclusions are as follows: 

• Despite stringent SI and SD limits, neutralino DM remains experimentally viable. The 
allowed regions of parameter space are large; they permit thermal, non-thermal, and multi- 
component neutralino DM, and include points with minimal fine-tuning. This result hinges 
on the proximity of viable parameter regions to numerous direct detection "blind spots" at 
which the SI or SD scattering cross-section vanish identically due to destructive interference 
in the couplings Ch^x or Czxx- have analytically identified all such cancellation points for 
neutralino DM and found that they require diverse sign choices among the mass parameters 
(summarized in Table 1). 

• Upcoming direct detection experiments will probe the bulk of neutralino DM parameter 
space, except for regions very close to the direct detection blind spots. Most of the regions 
left unscathed require fine-tuning between Mi, M2, and /x in order to track the blind 
spots. A notable exception to this is thermal bino/Higgsino and non-thermal Higgsino 
DM at low tan /3, which will evade both SI and SD experiments for the foreseeable future. 
Interestingly, this region of parameter space is theoretically favored: low tan /3 is required 
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by natural theories like ASUSY [10] in order to sufficiently boost the Higgs mass; it is also 
required by unnatural theories like Split Supersymmetry [9] , given Higgs mass constraints. 

• Bino/wino DM remains an attractive candidate for neutralino DM. As discussed in past 
works [11, 12, 13, 14, 15], the observed DM abundance can be accommodated by thermal 
freeze-out through coannihilations. We present here the ffist systematic analysis of the full 
parameter space of thermal, non-thermal, and multi-component bino/wino DM in relation 
to present and future experiments. 

Let us comment briefly on the relationship between our results and past work. Two recent 
studies of neutralino DM [16, 17] focus on the correlation between the size of the cross-section 
relevant for direct detection and naturalness in electroweak symmetry breaking. An earlier 
study [18] identified certain regions of parameter space with small direct detection cross-section. 
We note that these earlier papers all relied on scans over large parameter spaces, differing greatly 
from our approach of using simplified models. Well-tempered bino-Higgsino DM was studied in a 
simplified model approach in [19], and in the addendum it is claimed that the recent XENONIOO 
results exclude the case of thermal freeze-out. However, the case of a relative sign between fi 
and Ml, leading to destructive interference in the Higgs coupling Chxx^ i^o^ considered by 
these authors. 

The outline of our paper is as follows. In Sec. 2, we discuss the relevant experiments, focusing 
on present limits and future reach to probe the SI and SD cross-sections. We go on to discuss 
the cosmological history of the DM relic abundance in Sec. 3, including a review of the well- 
tempered neutralino. In Sec. 4, we identify regions of parameter space where the DM direct 
detection cross-section is suppressed. This suppression can come from purity of DM, which we 
discuss in Sec. 4.1, or from blind spots, where the Higgs or Z coupling vanish, which we classify 
in Sec. 4.2. We present detailed results on the present limits and reach for the simplified model 
where DM is a mixture of a bino and a Higgsino, in Sec. 5. Then we bring the wino into the 
spectrum in Sec. 6, and study DM that is a mixture of bino, wino, and Higgsino (with a special 
focus on the case where DM is dominantly bino/wino). Our conclusions appear in Sec. 7. 

We have included a few appendices with results that are of a more technical nature, but 
important for our study. Direct detection limits depend on the strange quark content of the 
nucleon, fg, and we use the lattice values [20] for our analysis. In App. A, we review the recent 
status of determinations of fg and quantify how our results are sensitive to this quantity. In 
App. B, we explain how we quantify the tuning of the DM abundance and direct detection cross- 
section, independently of a possible tuning of the electroweak scale. Throughout the paper, we 
use tree-level scattering cross-sections, and one may wonder how our limits, and specifically 
blind spots, are modified by loop corrections. In App. C we justify that these loop corrections 
should be small throughout most of the parameter space we consider. 

2 Observational Constraints 

In this section, we enumerate the experiments relevant to neutralino DM and broadly outline 
their status and future reach. However, before delving into the experimental limits we would 
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Figure 1: Present limits (filled or solid) and future reach (dashed) for SI/SD scattering of 
DM, shown in terms of the cross-section (left axis) or DM Higgs/Z coupling (right axis). For 
SI scattering we show the current limit from XENONIOO [1] as well as the projections for 
LUX [4], SuperCDMS [5], and XENONIT [3]. For SD scattering we show the current hmit 
from XENONIOO [6] on DM-neutron scattering, as well as the current limit from IceCube [2] 
on DM-proton scattering, assuming annihilations into W^W~ or ti (estimated). We also show 
our estimate for the reach of XENONIT [7] for DM-neutron scattering. 

like to ask: what is the characteristic size for the SI and SD cross-sections expected of neutralino 
DM which couples through the Higgs and Z bosons? Given the interactions, 

^ 3 ^ Kxx + x\^) + cz^x x^<y'xz,, (2) 

then in the limit in which the DM is heavier than the nucleon, the SI and SD cross-sections are 

asi = 8 X 10-^^ cm^ ' ago = 3 x lO'^^ cm^ ' . (3) 

While (TsD is typically considerably larger than ctsi, SI experimental constraints are commensu- 
rately stronger than SD, so these two limits are comparable in strength [21, 22]. Note that usi 
depends on nuclear form factors, in particular the strange quark content of the nucleon. For our 
analysis we adopt the lattice values of [20]. A more technical discussion of the strange quark 
content of the nucleon is contained in App. A. 

The SI scattering of DM with nucleons is highly constrained by null results from direct 
detection experiments. At the forefront of this experimental effort is XENONIOO [1], an un- 
derground, two-phase DM detection experiment which employs a 62 kg radio-pure liquid Xe 
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target. As shown to the left of Fig. (1), XENONIOO provides the current leading experimental 
limit on SI scattering; their latest limit uses an exposure of 0.02 ton x years. Throughout 
this paper, we present 90% CL. limits and reach, and we take the local DM density to be 
Po = 0.3 GeV/cm^, which is supported by a recent direct measurement using stellar kinematics, 
Po = 0.3 ±0.1 GeV/cm^ [23]. We do not attempt to incorporate astrophysical uncertainties into 
our analysis. 

XENONIOO will be far surpassed by XENONIT [3], which is projected to begin collecting 
data in ~ 2015 and should offer substantially better sensitivity due to a larger Xe target mass 
of ~ 2.2 tons. The projection in Fig. (1) shows the estimated limit with an exposure of 2.2 
ton X years. Meanwhile, LUX [4] — a direct detection experiment of similar design but with 
a Xe target volume of 350 kg — is slated for operation in ~ 2013, and has a projected reach 
between that of XENONIOO and XENONIT. In Fig. (1) we show the conservative LUX reach 
estimate from Ref. [4], which assumes an exposure of 0.08 ton x years and a light collection 
efficiency of 15%. A more optimistic light collection efficiency of 20% improves the limit by a 
factor of ~ 1.6. We also show the projected limit from SuperCDMS at SNOLAB [5], which 
is based on a complementary technology utilizing germanium detectors. Throughout the rest 
of the paper, we will focus on the XENONIOO limit, the conservative LUX estimate, and the 
XENONIT reach, but it is understood that the SuperCDMS reach would fall between the LUX 
and XENONIT estimates. We note that ton-scale Xenon detectors are also being pursued by 
LZ [24] and PandaX [25] and ton-scale Argon detectors are being pursued by the DarkSide 
collaboration [26]; we do not show their reach but our XENONIT curves should be taken as 
representative of the expected sensitivity of ton-scale liquid noble gas detectors. 

The SD scattering of DM with nucleons is constrained by direct detection experiments. The 
right of Fig. (1) shows the preliminary limit from XENONIOO, with 100 livedays, on the DM- 
neutron scattering cross-section [6] . We assume a nuclear shell model that leads to a conservative 
limit; a different model improves the limit by a factor of ~ 1.9, and this difference can be taken 
as an estimate of the uncertainty on the limit from nuclear physics. We also show an estimate of 
the reach with XENONIT. No official SD reach estimate has been released by the XENONIT 
collaboration, so we estimate the reach by rescaling the XENONIOO limit by the expected 
difference of exposure between the XENONIOO hmit (100.9 days x 48 kg) and XENONIT (2.2 
ton X years) [7]. 

There is also a constraint on SD DM-proton scattering from neutrino telescopes, which probe 
the annihilation into neutrinos of DM captured in the sun. We show the preliminary limit from 
IceCube [2] in Fig. (1), utilizing 79 strings and 317 days livetime. The limit from IceCube is a 
function of the neutrino spectrum, which depends on the DM annihilation products. We show 
the limit from DM annihilations into which is released by the IceCube collaboration, 

as well as our estimate of the limit if DM annihilates entirely into tt. In order to perform 
this estimate, we use the IceCube W~W~^ limit at fixed DM mass to determine the muon flux 
limit for mono-energetic W^s [27]. We determine the W energy spectrum resulting from top 
decays from annihilations at each DM mass using MadGraph [28], and estimate the upper limit 
by approximating the upper limit on the number of observed muons to be independent of the 
W energy. 
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We briefly comment on a few relevant constraints on neutralino DM other than direct detec- 
tion. There are indirect limits on DM annihilations into gamma rays. The strongest constraint 
comes from a combined Fermi analysis of 10 satellite galaxies using 2 years of data [29]. DM 
annihilating into W~W~^ is constrained to have a cross-section smaller than (av) < 10~^^ cm^/s, 
which as we will see places important constraints on DM with a non-thermal cosmology. The 
limit includes uncertainties on the haloes of the satellite galaxies, and should be viewed as 
conservative with regard to these uncertainties. In principle anti-proton measurements from 
PAMELA may set complementary limits [30]; however, we restrict our focus to Fermi because 
gamma rays, unlike antiprotons, are not sensitive to propagation uncertainties. There are also a 
few relevant limits coming from colliders. LEP2 constrains the charged components of Higgsinos 
and winos: m^± > 100 GeV [31]. There are now limits from the LHC constraining winos lighter 
than ~ 300 GeV that decay to a neutralino lighter than about 100 GeV [32, 33]. We will not 
consider the LHC limits further in this paper because we focus on DM heavier than 100 GeV, 
where these limits are not relevant. 



3 Relic Abundances and Well- Tempering 

WMAP observations are consistent with a relic abundance of DM given by [34]: 

^]^j^^/i2 = 0.111 ±0.006 (la). (4) 

Throughout our analysis, denotes the total relic abundance of neutralino DM, while Q^^'^ 
denotes the relic abundance of neutralino DM expected from thermal freeze-out alone. To be 
comprehensive, our analysis accommodates three scenarios for the cosmological history: 

• Thermal {Qohs = ~ ^x^^)- is solely comprised of neutralinos arising from thermal 
freeze-out. 

• Non-Thermal [Qohs = ^x^^)- DM is solely comprised of neutralinos, but thermal 
freeze-out either over- or under-produces. We assume that non-thermal processes either 
deplete or enhance the abundance to exactly saturate the WMAP constraint. 

• Multi-Component (i^obs > ~ ^x^^)- is partly comprised of neutralinos arising 
from thermal freeze-out. We assume that the balance of DM is provided by a secondary 
DM particle, e.g. axions. 

In the second and third cases, ^x^^ is not stringently constrained by WMAP measurements, 
so these scenarios offer greater freedom for evading experimental constraints. In the first case, 
however, the relic abundance is fixed to the observed WMAP value, and for x DM this typically 
requires a modest fine-tuning among parameters. This occurs because pure bino DM is over- 
abundant, while pure wino or Higgsino DM is under-abundant for masses below 1 TeV [35, 
36] and 2.7 TeV [36, 37], respectively. Thus, only a precise admixture of bino and wino or 
Higgsino — i.e. a well-tempered neutralino — can accommodate f^obs = ^x^^ [15] (for earlier refs., 
see [38, 39, 40]). 
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Figure 2: The thermal freeze-out abundance for x ~ h) (left) and x ~ {b, w) (right) is 
especially sensitive to parameters near the well-tempered cross-over region. The relic abun- 
dance is exponentially sensitive in bino/wino DM, where thermal freeze-out follows mostly from 
coannihilation. 



Fig. (2) shows the dependence of Vt^^'^ on Mi for the cases of x ~ (&, h) and x ~ (pi w). Here 
and throughout the paper, we compute relic densities with MicrOMEGAs 2.4.5 [41]. We do not 
include the effect of Sommerfeld enhancement, a non-perturbative effect which can substantially 
boost the annihilation cross-section of DM if it is much heavier than a force carrying particle. 
Sommerfeld enhancement through electroweak bosons is an especially important effect for wino 
DM > 2 TeV, which is not our focus. The parameters in Fig. (2) have been chosen so that 
in the limit of heavy Mi, x is dominantly a Higgsino or wino with mass 500 GeV. At low 
Ml, X is dominantly bino, and as Mi increases it gradually acquires a larger component of 
Higgsino/wino allowing it to annihilate more rapidly to final states involving W, Z, and h. 
As Ml approaches 500 GeV the abundance changes rapidly — partly because the mixing angle 
changes rapidly and partly because the LSP mass is approaching the mass of the next lightest 
neutralino and chargino states, allowing for coannihilation [42]. As is well-known, coannihilation 
is exponentially sensitive the mass difference between the DM and its neighboring states. The 
effective freeze-out cross-section for i = 1, 2, . . . , states is given by 

EN -n 

lfXi;;coann = T^^ 73 (5) 

= ( ^1 e-"("«/"^--i), (6) 



where x = fri-^/T and m-^^ = m^. Coannihilation dominates the transition region in the 
bino/wino case, leading to a curve that is much steeper than for the bino/Higgsino case. This 
is the consequence of the exponential sensitivity to mass parameters in Eq. (6). 

Because thermal neutralino DM requires a special relation among parameters, ~ |Mi| 
or IM2I, as we see in Fig. (2), it is worthwhile to pause and consider our philosophy on fine- 
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tuning. In this paper, we are agnostic about the possibihty of fine-tunings in both electroweak 

symmetry breaking and quantities relevant for DM phenomenology, such as and o"si,sd- 

We survey the entire parameter space of thermal and non-thermal neutralino DM, regardless of 
tuning. Indeed, it may be reasonable for a tuning of parameters to produce the observed fl if 
environmental selection plays a role in the DM abundance (although a tuning that produces a 
small a would be more surprising). In addition, we see in Fig. (2) that the relic density curves 
are steep for a wide range of Mi, as is expected from the phenomenon of well-tempering. In 
such a situation, where a large fraction of parameter space is highly sensitive to parameters, 
perhaps one should not be surprised to end up with parameters that appear tuned. Despite 
these misgivings about avoiding tuned regions, we do view it as interesting to identify regions of 
parameter space especially sensitive to parameter choice. At times we will quantify the amount 
of tuning in Q,^^^ and crgi go. In order to identify DM tunings independently of a possible 
electroweak tuning, we use a measure, defined in App. B, that mods out the dependence on the 
direction of parameter space that changes the electroweak VEV. 

4 Suppression of Dark Matter Scattering 

In general, SI scattering is mediated by squarks, Z bosons, or Higgses. Since squark-mediated 
scattering is model dependent — its effects become negligible for sufficiently heavy squark masses — 
we postpone our discussion of this scenario to Sec. 5.4. Similarly, we neglect scattering mediated 
by the heavy Higgses, which decouple quickly in the limit tua mz- This leaves scattering 
mediated by the light Higgs or Z, which may be suppressed compared to naive expectations by 
two effects. First, a suppression results whenever the DM is close to a pure gaugino or Higgsino, 
and second, the relevant amplitude exactly vanishes at critical values of parameters, which we 
call blind spots. 

4.1 Suppression from Purity 

The leading SI scattering is mediated by Higgs exchange, and the relevant coupling, Chxx^ orig- 
inates, at tree-level, from gaugino Yukawa couplings of the form h'^hb and h^hw. Hence SI 
scattering is suppressed if x is dominantly Higgsino or dominantly gaugino. Similarly, SD scat- 
tering does not occur for pure bino or pure wino because they do not couple to the Z, and 
likewise for pure Higgsino because as a Dirac fermion it carries no chiral couplings to the Z. 
Recall that the neutralino mass matrix takes the form. 



Since we are interested in Mi, M2, /i > Mz, Eq. (7) shows that gaugino- Higgsino mixing is gener- 
ically small, so that some cross-section suppression is expected for a typical point in parameter 



M^ = 




- ^g'v cos (3 \g'v sin /3 \ 

\gvcosj3 —\gvsm.j3 

-n 

-li 0. / 



(7) 



— \g'v cos (3 ^gv cos f3 
\ ^g'vsm(3 —^g'v cos (3 
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space. However, if the two lightest states are nearly degenerate the mixing between them, 6, 
can be appreciable, giving Chxx^ c^xx ^ with [15] , 

(sin /3 ± cos /3) sin 6^1^ / Mz 

for gaugino/Higgsino DM and 

^ _ sin2/3sin2^vK f M| 
~ 2 \fi{M2 - Ml) 

for bino/wino DM. Both results are valid for a mass splitting AM > Mz', 6w is the weak mixing 
angle and the signs in Eq. (8) refer to the cases fi ~ ±Mj. 

For successful thermal freeze-out with Q,^^^ = fiobs some degree of degeneracy is required, 
as seen in Fig. (2), so that SI and SD scattering may not be suppressed. However, significant 
suppression is expected for typical parameters in the cases of non-thermal or multi-component 
DM. 

So far we have been considering tree-level scattering, which vanishes for pure gaugino or 
Higgsino. But scattering between a pure Higgsino or wino and nuclei is generated by loop 
diagrams, for example 1-loop box diagrams with the exchange of two gauge bosons. Naively the 
1-loop scattering has a SI cross-section of a ~ IQ~(^'^-^^) cm^, which could be probed by the next 
generation of direct detection experiments. However, an accidental cancellation [43, 44] among 
various 1 and 2 loop diagrams leads, for pure Higgsino or wino, to cross-sections too small to 
probe at XENONIT, a < 10"^^ cm^. 





4.2 Suppression from Blind Spots 

To obtain a formally vanishing tree-level cross-section through purity, the gauginos or Higgsinos 
must be completely decoupled, Mi^2 or /i — >■ oo. We now consider a different possibility: special 
choices of parameters where the tree-level cross-section vanishes identically. At these blind spots, 
the gaugino and Higgsino masses are finite and the mixing is non-zero. 

Throughout our analysis, we assume that Mi, M2, and /i are real parameters, but carry 
arbitrary signs. However, only two of the three apparent signs are physical, as is clear from the 
field redefinition 

b ^ ib (10) 
w iw (11) 
hu,d -> -ihu,d, (12) 

which is equivalent to simultaneously sending the all the mass parameters Mi, M2, and fi to 
minus themselves. In many of our results, we will eliminate the unphysical, overall sign by fixing 
the sign of a single theory parameter to be positive. 

Let us denote the mass eigenvalues of M^ by m^.{v), where i = 1,2,3,4 and m^^ = is 
the DM mass. Here we have emphasized the explicit v dependence in the masses. The coupling 
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condition 


signs 


Ml 


Ml +/isin2/3 = 


sign(Mi//i) = -1 


M2 


M2 + /isin2/3 = 


sign(M2//i) = -1 


-fi 


tan /3 = 1 


sign(Mi,2//i) = -1* 


M2 


Ml = M2 


sign(Mi,2//i) = -1 



Table 1: Table of SI blind spots, which occur when the DM coupling to the Higgs vanishes 
at tree-level. The first and second columns indicate the DM mass and blind spot condition, 
respectively. All blind spots require relative signs among parameters, as emphasized in the 
third column. *For the third row, the blind spot requires that /i and Mi (M2) have opposite 
signs when M2 (Mi) is heavy. 

of any of neutralino to the Higgs boson can then be obtained by replacing v v + h, as dictated 
by low-energy Higgs theorems [45, 46]: 

Axx = \^x^i^ + h)XiXi (13) 

= ^^x.(^)xa. + ^^^^te. + C?(/^'), (14) 

which implies that dm^^{v)/dv = Ch^iXi [47, 48]. 

Consider the characteristic equation satisfied by one of the eigenvalues rriy^-iy), 

dei{M^ - lm^^{v)) = 0. (15) 

Differentiating the left-hand side with respect to v and setting dm^.{v)/dv = ChxiXi = 

then obtains a new equation which defines when the neutralino of mass ?^xi(f) ^ vanishing 

coupling to the Higgs boson^: 

(m^,(t;)+/xsin2/3) (^mx^(t;)-^(Mi + M2 + cos2^v^(Mi-M2))^ = 0. (16) 

The above equation implies that for regions in which Ch^iXi ~ ^1 ^xX'^) entirely independent 
of V. At such cancellation points, mjx-(f) = m^. (0), so the neutralino mass is equal to the mass 
of a pure gaugino or Higgsino state and m^.(t>) = Mi, M2, — /i. As long as Eq. (16) holds for the 
LSP mass, m^^(f), then the DM will have a vanishing coupling to the Higgs boson, yielding a 
SI scattering blind spot. It is a nontrivial condition that Eq. (16) holds for the LSP, rather than 
a heavier neutralino, because for some choices of parameters the DM retains a coupling to the 
Higgs but one of the heavier neutralinos does not. We have identified these physically irrelevant 
points and eliminated them from consideration. The remaining points are the SI scattering 

^We have checked that Eq. 16 can also be derived using analytical expressions for bilinears of the neutralino 
diagonalization matrix from Ref. [49]. 
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blind spots for neutralino DM, 



spin-independent 
blind spots 



m 



m 



Xi 



Ml, M2, —fjL, and m^^ + /isin2/3 = 
Ml = M2, 



(17) 



where in the first line, m^^ = Mi,M2, — /i, depending on whether the LSP becomes pure bino, 
wino, or Higgsino, respectively, in the f — )■ limit. Note that the blind spots in Eq. (17) only 
appear for certain choices of relative signs. In the first line, for example, if m^^ = Mi(M2), then 
H and Ml (M2) must have opposite signs; when m^^ = —fi, then /i must have the opposite sign 
of Mi(M2) when M2(Mi) is heavy. For the second line, the blind spot occurs if fi and Mi = M2 
have opposite signs. The complete set of conditions required for a SI blind spot are summarized 
in Table 1. 

Destructive interference between light and heavy Higgs exchange may also produce cancel- 
lations in the SI cross-section [50], but these are outside the scope of this work. We consider 
interference between Higgs and squark exchange in section 5.4. 

Next, let us consider SD scattering, which is mediated by Z boson exchange. The coefficient 
of the relevant operator vanishes for neutralino DM when 



yielding a blind spot for SD direct detection. The cancellation of the SD Z boson coupling to 
DM can be understood from symmetry arguments: when Vu = Vd, the DM Lagrangian enjoys 
an enhanced symmetry under which u d. In this limit left-right parity is restored and hence 
the parity-violating Z coupling which mediates SD scattering will vanish. 

So far our discussion of blind spots has been tree-level. One may wonder how the blind spots 
change when loop corrections are included. Loop corrections have not been computed in the 
full parameter space, but only for the simplifying assumption of pure DM [43, 44], as discussed 
above. But our expectation is that the loop corrections are small, generically resulting in a small 
shift in the location of the blind spots. Moreover, at a typical point in parameter space, the 
mixing angles are small and the multiloop result for pure Higgsino or wino will approximately 
apply, leading to a cross-section too small to probe in upcoming experiments like XENONIT. 
Full consideration of loop corrections is beyond the scope of our study, but we estimate the size 
of these corrections in App. C. 

5 Bino/Higgsino Dark Matter 

In this section we consider the present and future status of non-thermal, multi-component and 
thermal bino/Higgsino DM. Mixed bino/Higgsino has been studied in a variety of contexts 
and more recently has been re-examined in light of results from direct detection experiments 
[16, 17, 19]. Here we take a simplified model approach to bino/Higgsino DM, decoupling all 
superpartners, other than the bino and Higgsinos, and all Higgs-like scalars other than the SM- 
like state near 125 GeV. Thus DM is described by just three parameters, (Mi, /x, tan/3). Our 



spin-dependent 
blind spot 



tan/3 = 1, 



(18) 
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analysis applies to this decoupled limit of the MSSM, NMSSM and to Split Supersymmetry. 
At the end of the section we consider additional effects that arise when the squarks are not 
decoupled. Some effects from non-decoupling of the wino are illustrated in the next section. 
For simplicity we remove a physical phase by imposing CP conservation on the neutralino mass 
parameters, but we study the effect of the physical sign between /i and Mi. Our convention is 
to take P in the first quadrant and choose Mi positive, allowing both signs of fi. Our numerical 
results here, and in Sec. 6, use MicrOMEGAs 2.45 for cross-sections [51] and relic densities [41]. 

We may understand the results of the following subsections by considering the variation of 
the SI and SD elastic scattering cross-sections within the bino/Higgsino parameter space, as 
shown in Fig. (3), arising from the tree-level exchange of a SM Higgs of mass 125 GeV for SI 
and from the Z boson for SD. At large tan /3 (lower panel) the SI and SD cross-sections become 
independent of the sign of fi. This happens because in this regime the bino mixes negligibly 
with the down-type Higgsino and so the sign of fi can be removed by a field redefinition and 
is unphysical. Both the SI and SD contours fall off with increasing /i or Mi when the other 
parameter is kept fixed, as expected from the vanishing of the bino/Higgsino mixing angle. 
As fi approaches Mi, the mixing angle maximizes, and the contours show a ridge along the 
line Ml ~ /i, with the cross-sections dropping off steeply on both sides. This behavior can be 
understood from the discussion in section 4.1; the ridge corresponds to the region with large 
mixing between the lightest two nearly degenerate states, given by eq. (8), from which it follows 
that the ridge becomes steeper at large masses. Furthermore, this region of maximal scattering 
cross-section coincides with the well-tempered line, since large mixing is also necessary to achieve 
the observed relic abundance. 

The black dashed lines show blind spots for SI scattering with Chxx = 0? arising from the 
relation Mi -|- sin 2/3 /i = 0. For high tan/3 this occurs in a region where the SI cross-section 
is highly suppressed by a small mixing angle, but at low tan (3 (upper panel) the blind spot 
cuts a gorge in the ridge at negative /i. As we will see, the proximity of this blind spot to the 
region with large mixing angles has important implications for the observability of thermal DM, 
although the rapid variation of the contours implies an enhanced tuning of the cross-section. 



5.1 Non-thermal Dark Matter with = Qobs 

We begin by considering the limits on bino/Higgsino DM without imposing that thermal freeze- 
out provide the observed relic abundance. Fig. (4) depicts contours of il^^V^obs — the ratio 
of the thermal yield of neutralino DM to the observed relic abundance — together with current 
experimental constraints in the (/i. Mi) plane at tan/3 = 2,20. For |Mi| ^ DM is bino-like 
and flx^^ is over- abundant; for |Mi| ^ DM is Higgsino-like and flx^^ is under- abundant. In 
these regimes we have evaluated constraints assuming a non-standard cosmological history in 
which entropy production or non-thermal DM production, respectively, ensures a final neutralino 
abundance of f2obs = ^ ^x^^- 

According to Fig. (4), thermal bino/Higgsino DM at low tan/3 is excluded up to rrix — 800 
GeV for /i > but practically unconstrained for fi < 0. At high tan P, however, thermal 
bino/Higgsino DM is excluded for m-^ ~ 500 GeV for either sign of /i. Meanwhile, non-thermal 



14 



SI and SD cross-sections for b/h 

tan yS = 2 




yu [GeV] 



tan yS = 20 




yu [GeV] 



Figure 3: Contours of the tree-level cross-sections for SI (solid red) and SD (dashed blue) 
scattering of bino/Higgsino DM. The brown band denotes regions with ^x^^ within ±3cr of ^ohs- 
The black dashed line is the SI blind spot, Ch^x = 0, arising from the relation Mi + /i sin 2/3 = 0. 
The central gray region is excluded by LEP. 

bino-like or Higgsino-like DM is, at present, rather poorly constrained by direct detection on 
account of the relatively small mixing, and therefore small couplings to the Higgs and Z. Con- 
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Figure 4: Current limits on bino/Higgsino DM with = Qohs for tan/3 = 2 (upper), 20 
(lower). Dotted brown lines are contours of Q^^^ /Qohs, and the brown band shows the region 
having ^Ix^^ within ±3cr of fiobs- Regions above (below) the brown band require an enhancement 
(dilution) of the DM abundance after freeze-out. Regions currently excluded by XENONIOO, 
IceCube, Fermi, and LEP are shaded. The black dashed line is the SI blind spot, Ch^x = 0, and 
is close to (far from) the brown band for low (high) tan/3. 
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versely, indirect detection limits do exclude non-thermal Higgsino-like DM above the W^W~ 
threshold up to yU ~ 250 GeV, but well-tempered DM, which has a smaller annihilation cross- 
section, evades this bound. 

The Fermi limit on DM annihilation to W^W~ primarily comes from constraints on photons 
created in the decays of hadrons. In order to obtain the Fermi exclusion region shown in Fig. (4), 
we include the annihilation cross-sections to both W^W~ and ZZ, weighted by the relative 
hadronic branching ratios. 

The structure of Fig. (4) follows from the fact that the leading experimental constraint on 
bino/Higgsino DM is on SI scattering at XENONIOO. In general, thermal neutralino DM tends 
to be the most constrained by SI direct detection, simply because DM carries an 0{1) fraction of 
bino and Higgsino that furnishes a non- vanishing coupling to the Higgs. In contrast, parameter 
regions corresponding to non-thermal DM are more weakly constrained, since pure bino and pure 
Higgsino DM do not couple directly the Higgs boson. That said, even well mixed neutralino DM 
can be decoupled from the Higgs if the theory parameters reside on the SI blind spot defined in 
Eq. (17), 

Ch-^^ oc Ml + /isin2/3 = 0, (19) 

which is only allowed for /i < 0. At low tan /3, the SI blind spot occurs near the well-tempering 
region, \Mi\ ^ The tan/3 = 20 plot in Fig. (4) is approximately symmetric under O — /x, 
for the reasons noted earlier. 

At present, limits on the SD scattering cross-section are dominated by IceCube bounds. 
These provide a complementary constraint for lighter DM near the SI blind spots, since SD 
scattering cross-sections are unaffected by the vanishing of the coupling to the Higgs. As dis- 
cussed in Sec. 2, IceCube provides bounds on DM annihilation to W~^W~] in order to generate 
the exclusion regions in Fig. (4) and Fig. (5), we compare these bounds to the DM annihilation 
cross-section into W^W~, Zh, and ZZ, weighted by their branching ratio to neutrinos relative 
to W^W~ . Previous studies [52, 53] have also considered limits on neutralino DM annihilation 
from IceCube, using older data or projections. 

Fig. (5) is identical to Fig. (4) except it depicts projected reach instead of current limits. 
Comparing Fig. (5) and Fig. (4), LUX and XENONIT will provide a very powerful probe of 
both thermal and non-thermal bino/Higgsino DM. Currently only narrow wedges of the (/x. Mi) 
plane are excluded. These wedges lie along the thermal band, but even the exclusion of some 
thermal regions is marginal and subject to astrophysical uncertainties. Over the next few years, 
LUX and XENONIT will explore most of the parameter space with DM masses up to 1 TeV, 
and much of the region up to 2 TeV, offering a remarkable opportunity for discovery. If no signal 
is seen, LUX will exclude a large fraction of thermal bino/Higgsino DM and, XENONIT will 
exclude the entire parameter space of thermal bino/Higgsino DM for tan/3 > 2, except for the 
case of almost pure Higgsino. Interesting blind spot regions remain for lower tan /3, as discussed 
in the following subsections. 

On the other hand, even in the absence of a signal, significant parameter regions for the non- 
thermal case will remain. Bino-like DM is permitted for /i < near the SI blind spot for bino-like 
DM defined in Eq. (19). Meanwhile, non-thermal Higgsino-like DM is highly unconstrained at 
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Figure 5: Same as Fig. (4) except for future reach rather than current hmits. The dashed green 
hues show the projected SI reach of LUX, while the shaded regions give the projected reach 
for XENONIT, both SI and SD. The shaded cyan region is the current Fermi exclusion, as in 
Fig. (4). 
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Figure 6: Limits and projected reaches for multi-component bino/Higgsino DM with fl^ = fl\^ . 
Dotted brown hues are contours of 0,^^^ /Q,ohs for tan/3 = 2. The light gray regions are excluded 
by overabundance of neutralino DM, while the edge of this region has ^x^"^ = ^ohs- In the 
remainder of the plane x is just one component of multi-component DM. The present limit 
from XENONIOO is shown shaded, while the projected reaches of LUX and XENONIT, both 
SI and SD, are shown as dashed lines. 

low tan (3 because it corresponds to the SI blind spot for Higgsino-like DM in Eq. (17), 

Chxx oc -1 + sin 2/3 = 0. (20) 

Some of these allowed regions will be probed by experiments sensitive to the SD scattering cross- 
section. Intriguingly, the case of non-thermal Higgsino DM at low tan /3 resides simultaneously 
in a blind spot for SI and SD scattering! Furthermore, this region allows low values of fi, 
and therefore relatively natural theories of electroweak symmetry breaking. In addition, large 
unnatural regions with /i > 1 - 2 TeV will remain viable, but require late entropy production, 
especially for low Mi. 

5.2 Multi-Component Dark Matter with = Q^^^^ < ilobs 

Here we repeat the analysis of the previous section under the assumption that the present day 
relic abundance of neutralino DM is given by = ^x^\ with the balance of cosmological 
DM arising from some other source. Fig. (6) depicts both the current limits and the projected 
reach for such multi-component neutralino DM, for tan/3 = 2. Region shaded light gray have 
Qx^^ /Qohs > 1 and are thus excluded, while regions with Qx^'^ /Qohs < 1 have a depleted abun- 
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dance of neutralino DM. Direct detection limits are then ameliorated, since the rate of WIMP- 
nucleon scattering is proportional to the incident flux of DM particles, and thus to Q^^^ /Qohs- 
DM annihilation cross-sections are suppressed by the square of this ratio, making indirect de- 
tection limits irrelevant. The edge of the light gray shaded region has flx^^ = ^ohs and therefore 
has thermal bino/Higgsino DM, which is explored for all tan/3 in the next sub-section. 

It would appear to require a coincidence for otherwise unrelated stable particles to have 
comparable relic abundances, as would be needed for multi-component neutralino DM with an 
0{1) fraction of f2obs- However, environmental selection may provide a possible explanation 
[54, 55]. While the expected abundances in the various components depend on multiverse 
distribution functions, it is likely that they are very roughly comparable. Alternatively, multiple 
sectors of a theory may participate in the WIMP miracle, in which case each WIMP would 
independently attain a thermal abundance near Qohs- Consider for example the case that the 
bino/Higgsino makes up a third of the total DM. Current limits from XENONIOO and IceCube 
allow any mass above the LEP exclusion limit. However, XENONIT will probe a large and 
interesting region, pushing the mass up to about 500 GeV if no signal is seen. Much of the 
remaining space is dominantly Higgsino. 

5.3 Thermal Dark Matter with = Q^^^^ = Hobs 

We now restrict our analysis to well-tempered bino/Higgsino DM. We fix the value of Mi using 
the relic abundance constraint, flx^"^ = ^ohs, reducing the parameter space by one dimension. 
Thus we can show the entire parameter space relevant for thermal bino/Higgsino DM in the 
(/i, tan P) plane. Fig. (7) depicts current limits and projected reach for thermal neutralino DM 
as a function of these parameters. For /i > 0, SI direct detection currently rules out thermal 
bino/Higgsino DM for /i < 650 — 800 GeV, depending on tan/3, fi < 0, however, is almost 
completely unconstrained, except for a small region around the tt threshold for DM annihilation 
at high tan /3. 

Future direct detection experiments will cover the entire well-tempered parameter space for 
fi > 0, and almost all of it for /i < 0, with the exception of a region around the blind spot 
cancellation given by Eq. (19). The DM coupling to the Z does not vanish at the SI blind 
spot, however, so that SD direct detection will set complementary limits in this region, with 
XENONIT probing up to /i ~ -500 GeV. 

Fig. (8) quantifies the fine-tuning of thermal bino/Higgsino DM. These plots depict the 
sensitivity of the relic density and the SI scattering cross-section with respect to the ultraviolet 
parameters in the theory — namely. Mi, M2, /i, m^^, m^^, 5/i defined at the weak scale — using 

the measure defined in App. B. As might be expected from the steepness of the flx^^ curve in 
Fig. (2), well-tempering requires (9(10%) tuning throughout the parameter space. A similar, 
relatively mild tuning of the direct detection cross-section is required to evade the upcoming 
LUX bounds on SI scattering. Evading the limits from XENONIT, however, requires a cross- 
section tuning of at least 1%, suggesting that it would be unlikely for thermal bino/Higgsino 
DM to remain hidden. 

This tuning of the cross-section, however, is directly correlated with that required to obtain 
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Figure 7: Current limits from XENONIOO and IceCube (top), expected reach of LUX and 
IceCube (middle), and expected reach of XENONIT (bottom) for SI and SD scattering, fixing 
Ml at every point to accommodate fl^^'^ = Qohs- The black dashed line is the blindspot, Chxx = 0- 
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Figure 8: Relic density and SI cross-section tuning for well-tempered bino/Higgsino, with the 
reach of LUX and XENONIT also shown in the second panel. Tuning of the relic density is 
typically between 2 - 10%. If XENONIT does not see a signal, tuning of the SI cross-section 
will be < 1%. The interpretation of this as being unnatural is unclear however, as most of the 
region of 1 < tan/3 < 2 with < has large A„. We describe our methodology for computing 
tuning in App. B. 

the correct relic abundance at low tan/3, since both a large mixing angle and a small Higgs 
coupling require \fi\ ~ \Mi\. Furthermore, many theories, both natural and unnatural, require 
small tan /3 in order to explain the 125 GeV Higgs mass, as in ASUSY and Split Supersymmetry. 
Thus the region of parameter space which evades XENONIT is exactly the same region in which 
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one might expect to find oneself given both the Higgs mass and the observed rehc abundance. 
In this case, perhaps a large tuning of the cross-section should be unsurprising, since such a 
tuning is generic within the allowed parameter space. 

5.4 Squark Effects 

Our analysis thus far has neglected the effects of squarks. For the natural theories, however, 
it is reasonable to consider squark masses and /i that are not exceedingly large, so as not to 
exacerbate fine-tuning in electroweak symmetry breaking. In many such scenarios, the DM 
is mixed bino/Higgsino, and the effects of light squarks can play an important role on the 
physics [56, 57, 58]. 

Fig. (9) depicts present limits and future reach for bino/Higgsino DM including the effects of 
light squarks at tan/3 = 20 and /i < 0. The analogous constraints for fi > are more stringent. 
Here we have chosen degenerate 1st and 2nd generation squarks at m?, with the 3rd generation 
decoupled for simplicity. Squarks of the first generation tend to have the biggest effect because 
their exchange allows the bino to couple directly to valence quarks. At each point in the plot 
we have fixed Mi to accommodate thermal DM, Q.^^^^ = fiobs- In the limit of m? — )■ oo, this plot 
asymptotes to the current XENONIOO limit shown for bino/Higgsino DM in Fig. (7). 

The bino/Higgsino /squark space of Fig. (9) can be divided into two regions, depending on 
whether the squarks are heavier, or lighter than the Higgsino. When the squarks are lighter than 
the Higgsino, rrig < the correct abundance follows from squark/bino coannihilation. We find 
that this entire region of Fig. (9) is ruled out by XENONIOO. When rriq > the abundance 
follows from bino/Higgsino well-tempering, as in Figs. 2 and 7. In this region, we see that 
as the squark masses drop the XENONIOO limit disappears. There are two important effects 
contributing to this, (1) the contribution to the SI scattering amplitude from s-channel squark 
exchange destructively interferes with the contribution from t-channel Higgs boson exchange, 
and (2) squark diagrams increase the annihilation rate, leading to a smaller bino-Higgsino mixing 
angle and a smaller DM-DM-Higgs coupling, when the relic density is fixed. 

Interestingly, as the squark mass approaches the LSP mass, limits from supersymmetry 
searches at the LHC are also alleviated. In particular, we have plotted the limit on the 
squark/bino simplified model of [59]. The purple dashed line represents the limit, presented 
by ATLAS, when the gluino is decoupled. The purple shaded region represents our estimate 
of the limit, when rUg = 2 TeV, which is stronger due enhanced squark production with a t- 
channel gluino. In order to estimate this limit, we have reproduced the ATLAS search using 
Pythia 6.4 [60] for event generation, PGS for crude detector simulation [61] and NLLfast for 
the NLO [62] and NLL [63] squark cross-sections. Both constraints are somewhat weakened 
near the region of bino/squark degeneracy. Regarding future reach. Fig. (9) indicates that the 
projected SI and SD constraints from XENON and LUX rule out a large fraction of the allowed 
space of bino/Higgsino DM with a light squark. 
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Figure 9: Impact of squarks on thermal bino/Higgsino DM, with < and tan /3 = 20. At each 
point Ml has been chosen so that VL^^^ = flohs, except in the gray region where freeze-out always 
yields overclosure. The upper left region, where freeze-out is dominated by squark-neutralino 
coannihilation, is excluded by XENONIOO. However, in the lower right region the XENONIOO 
limit becomes less powerful as the s-channel squark exchange amplitude has the opposite sign 
to the t-channel Higgs exchange diagram. The purple region is excluded by an LHC search for 
jets and missing transverse energy, with the gluino mass fixed at 2 TeV. This ATLAS search 
becomes less powerful as the gluino mass is increased, and the excluded region becomes bounded 
by the purple dashed line if the gluino is decoupled. The currently allowed region, shown in 
white, mostly has a SI scattering cross-section that is not far below the current bound, so that 
LUX will have a large discovery potential. In the absence of a signal at LUX (XENONIT) the 
only surviving region will be the narrow band between the dashed green (red) lines. 

6 Bino/Wino(/Higgsino) Dark Matter 

We now consider the effects of including the wino in the spectrum. Compared to the previous 
section, reintroducing the wino adds an extra parameter, so that now we have a four dimensional 
parameter space of (Mi, M2, /i, tan In general, the LSP is now a combination of bino, Higgsino 
and wino, but much of our attention will focus on the case of a dominant bino/wino mixture. 
Even when the dark matter has a very small Higgsino component, the value of the parameter is 
crucial for direct detection: in the limit of decoupled /i, bino/wino dark matter has vanishingly 
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Figure 10: Phase diagram of neutralino DM in the (Mi,M2) plane, keeping /x fixed and less 
than 1 TeV. The red, green, and blue regions correspond to DM that is mostly bino, wino, 
or Higgsino-like. The thermal abundance, equals the observed abundance, f^obs, along 

the brown curve, which resides at the boundary of the bino region and wino/Higgsino regions. 
Within the bino-like region, the thermal abundance is too large and dilution is required; within 
the wino and Higgsino regions the thermal abundance is too small and additional neutralino 
production is required. 

small couplings to the Higgs and Z bosons. While mixed bino/wino dark matter has been 
explored in the past [11, 12, 13, 14, 15, 64], it has received substantially less attention than other 
limits of neutralino dark matter. Part of the reason for this neglect is theoretical prejudice. In 
particular, since thermal bino/wino dark matter originates via coannihilation, working models 
typically require Mi ~ M2, which is disfavored by gaugino unification. Moreover, as discussed 
in Sec. 3, the coannihilation region is exponentially sensitive to the mass splittings in the theory. 
Obviously, non-thermal or multi-component bino/wino dark matter require no such constraint 
on the masses, and have greater freedom to evade bounds. 

In this section we present a detailed study of non-thermal, mult i- component, and thermal 
bino/wino DM, focusing on present limits and future reach. Once again, we remove the physical 
phases in the neutralino parameters by assuming CP conservation. In contrast with the previous 
section, however, there are now two physical, relative signs in the theory. We continue to take /3 
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in the first quadrant, and we fix // > for non-thermal and multi-component DM. For thermal 

bino/wino DM, however, we use the constraint on relic density to fix Mi, which we take to be 
positive, allowing /i and M2 to have either sign. 

Fig. (10) shows a schematic slice of the parameter space relevant for bino/wino DM at fixed 
and tan^. For IM12I < the dark matter is dominantly either bino-like or wino-like 
depending on the hierarchy between \M]\ and IM2I. Even along the well-tempered line, dark 
matter is mostly bino, since the relic abundance is set by coannihilation with small mixing 
angles. As |Mi| approaches at large IM2I, the bino(wino)/Higgsino mixing angle increases 
until we recover the well-tempered DM considered in the previous section. In the limit of 
large |Mi|, we may have a mixed wino/Higgsino LSP for IM2I ~ but the observed relic 
abundance cannot be achieved for IM2I < 1 TeV due to the large annihilation cross-section 
to W^W~ . The final possibility, also considered in previous sections, is a dominantly Higgsino 
LSP if < |Mi^2|, which has an over-abundant (under-abundant) thermal relic density for 
> 1 TeV (1/^1 < 1 TeV). In this section we will focus on the upper and lower triangles 
of Fig. (10), containing non-thermal or multi-component wino DM and bounded by thermal 
bino/wino DM. 

6.1 Non-thermal Dark Matter with Qy. = r^obs 

Fig. (11) shows the current limits and future reach for neutralino dark matter in a slice of 
the full parameter space, in particular in the (M2,Mi) plane at tan^ = 2 and // = 750 GeV. 
The narrow brown bands in these figures show regions with ^x^"* within 3a of f^obs- The four 
physically independent sign choices are realized by allowing Mi 2 to range over both positive 
and negative values, yielding four quadrants. The relic thermal abundance depends largely on 
gauge interactions, which are independent of the sign choice, so the brown bands are largely 
the same in the four quadrants. The horizontal brown bands occur at |Mi| = mh/2 correspond 
to annihilation of mainly bino dark matter through the Higgs pole. The brown bands with 
slope near 45° have |Mi| ~ IM2I < correspond to bino/wino dark matter, and are narrow 
because coannihilation is operative. As IM2I increase these bands fiatten out and represent 
bino/Higgsino dark matter as the wino decouples. Here the bands are thicker as coannihilation 
is no longer present. The nature of the LSP varies across the (M2, Mi) plane of Fig. (11) as 
described in the beginning of this section. At low IM2I there is a light chargino that is excluded 
by LEP, as shown by the light gray band. The dark gray area is also excluded because the LSP 
is a chargino. 

The dark matter direct detection limits and reaches of Fig. (11) assume that — ^obs 
throughout the plane. This requires that, for any M2, regions with |Mi| larger (smaller) than 
that giving the brown bino/wino or bino/Higgsino band requires enhancement (dilution) of the 
dark matter abundance after freeze-out. The shaded red region of Fig. (11) is excluded by the 
recent XENONIOO search for SI scattering of dark matter. This region is highly asymmetric in 
the four quadrants because direct detection depends on the Higgs coupling Chxx which depends 
critically on the sign choices. In particular, the blind spots for SI scattering of Eq. (17), with 
Chxx — 0' occur when Mi = —//sin 2^ < IM2I and M2 = —//sin 2^ < |Mi|, as shown by the 
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Figure 11: Current limits and future reach for dark matter with Vt^ = ^lohs- The brown 
band shows the region having ^Ix^^ within 3a of fiobs- Regions with \Mi\ larger (smaller) than 
the brown band require an enhancement (dilution) of the dark matter abundance after freeze- 
out. Regions currently excluded by XENONIOO and FERMI are shaded in the left-hand figure, 
and the projected SI reach is depicted with dashed green lines for LUX and in shaded red for 
XENONIT on the right. The horizontal (vertical) black dashed lines are the SI blind spot, 
Chxx = 0) for bino-like (wino-like) dark matter. 

horizontal and vertical dashed black lines in Fig. (11), respectively. A third blind spot occurs 
with > Ml = M2 > —/i sin 2/3. It is only the quadrant which has no blind spot (Mi, M2 > 0) 
that is substantially constrained by current direct detection bounds from XENONIOO, being 
largely excluded except for a small region around the tt threshold. On the other hand, the 
remaining sign combinations of Mi and M2 are more or less unconstrained. 

Thermal dark matter is excluded for > 500 GeV for Mi,M2 > 0. In contrast with the 
limits from bino/Higgsino dark matter, limits on bino/wino DM actually become more stringent 
at larger values of Mi and M2. This behavior arises simply because SI scattering from Higgs 
boson exchange for bino/wino DM is mediated by the Higgsino fraction of the LSP. As Mi 
and M2 increase, then at fixed fi, the dark matter acquires a larger Higgsino fraction and SI 
scattering is proportionally larger. Present limits on thermal bino/wino DM are weak because 
/X can be modestly decoupled while still accommodating a thermal relic abundance. This is 
possible because the dominant process during thermal freeze-out is coannihilation of bino-like 
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Figure 12: Limits and projected reaches for multi-component dark matter with Vt^ = 
for tan/3 = 2, /i = 750 GeV. Dotted brown lines are contours of Q^^^ /Qohs- The gray regions 
are excluded by overabundance of neutralino dark matter, while the edge of this region has 
^Ix^^ = f^obs- In the remainder of the plane x is just one component of multi-component dark 
matter. The present limit from XENONIOO is shown shaded, while the projected reaches of 
LUX and XENONIT, both SI and SD, are shown as dashed lines. 

dark matter with charged and neutral winos. Even for larger values of /i than shown in Fig. (11), 
the dark matter can be thermally equilibrated with the winos, permitting efficient coannihilation. 

Finally, non-thermal wino-like DM is also constrained by Fermi gamma ray searches. The 
blue region of Fig. (11) shows that dark matter being composed of winos is excluded up to a mass 
of about 500 GeV. Compared to the limits described in Sec. 5, indirect detection constraints are 
more stringent for wino-like DM than they were for Higgsinos due to group-theoretic factors. 
As before, Fermi is currently unable to constrain thermal dark matter, which has a smaller 
annihilation cross-section. For prior applications of Fermi limits to non-thermal wino DM, see 
for example [65, 66, 67]. 

The right-hand plot of Fig. (11) is the same as the left except it shows future reach rather 
than current limits. Here we see that XENONIT will exclude large swaths of dark matter, 
leaving small patches corresponding to SI blind spots (aside from a small region of very light 
non-thermal bino DM, which we will not consider further). In particular, the horizontal and 
vertical black dashed lines denote the blind spots for bino-like and wino-like DM, respectively. 
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Consistent with Eq. (17), for larger values of tan ^5, these bhnd spots arise at smaller and smaller 
values of \Mi\ and \M2\, respectively. Spin-dependent direct detection limits, while not absent, 
are subdominant to SI limits throughout the plane of Fig. (11) and are therefore not shown. 

The contrast between the right and left plots is stark: currently TeV-scale neutralino dark 
matter is poorly constrained by direct detection experiments; but over the coming few years we 
can expect a much deeper probe of the theory yielding a large discovery region. The absence 
of a signal will require either that parameters lie close to a blind spot or that mass parameters 
are above the TeV scale. These conclusions persist with slices through the parameter space at 
other values of /i and tan/3. 

The diagonal black dashed line denotes the Mi — M2 bhnd spot for bino/wino DM. In- 
terestingly, this cancellation region coincides with the parameter space which accommodates 
a thermal relic abundance. In the event that the bino/wino DM abundance is anthropically 
selected, the fine-tuning imposed to acquire the correct relic abundance automatically induces a 
cancellation in the SI scattering cross-section of dark matter — no additional tuning is required. 

6.2 Multi- Component Dark Matter with = Q^^^^ < ilobs 

Fig. (12) shows the same slice through parameter space as Figure (11), so that the nature of the 
LSP in the various regions is as before. The key difference is that we now assume that the relic 
LSP abundance is given purely by thermal freeze-out. Hence, in addition to regions excluded 
by LEP searches for charginos and chargino LSPs, there is also a large region excluded by the 
overproduction of dark matter — essentially the entirety of the space containing a dominantly 
bino LSP. One quadrant is almost entirely excluded by these considerations, although a narrow 
strip close to thermal dark matter is allowed. 

Contours of il^^^ /riobs are shown as dotted brown lines. Regions shaded red have been ex- 
cluded by XENONIOO, while future projected reaches are shown by dashed lines. Both the hmits 
and the reaches are highly asymmetrical in the four quadrants, which can be understood from 
the locations of the blind spots. For example, XENONIT can limit the fraction of neutralino 
dark matter to less than 1% for Mi^2 > 0, but does not reach 10% if M2 is negative. 

6.3 Thermal Dark Matter with = U^^^^ = Qohs 

We now restrict our attention to thermal bino/wino dark matter, which once again allows us to 
reduce the parameter space of interest by one dimension. As before, we solve for Mi using the 
thermal relic constant, assuming Mi > 0. We consider slices of the remaining parameter space 
at fixed fjL and fixed tan /3. Fig. (13) depicts limits and reach for thermal bino/wino dark matter 
in the (M2, tan/3) plane at /i — —750 GeV. The SI blind spot for bino- like dark matter is shown 
in the plot using the expression from Eq. (19). 

As discussed in the previous section, there is no current limit from XENONIOO given this 
choice of the sign of /x. Furthermore, both LUX SI and XENONIT SD limits will only constrain 
the regions of parameter space with IM2I ^ which have a mixed bino/Higgsino LSP. Spin- 
independent limits from XENONIT will cover much of the space with M2 < 0, with the exception 
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Figure 13: Limit and reach for thermal bino/wino in the M2,tan/3 plane for /i = —750 GeV. 
Ml is fixed by requiring Vl^^ = flohs- The black dotted lines correspond to the SI blind spot for 
bino-like DM given in Eq. (19). The XENONIT SI and SD exclusion reach is shown shaded in 
red and blue, respectively, while the LUX SI reach is shown with a dashed green line. The LEP 
chargino exclusion is shaded in gray. 

of a region around the SI blind spot. For M2 > 0, however, the reach is weakened by virtue of 
the proximity of the Mi = M2 blind spot. Note that no additional tuning is required beyond 
that which is needed to get the correct thermal relic abundance, only a discrete choice of sign. 

Finally, Fig. (14) shows the current limit and expected reach for thermal bino/wino dark 
matter in the (M2, /x) plane at tan /3 = 2. Because of the location of the SI blind spots. Fig. (14) 
depicts much weaker constraints for negative fi than positive /i. For small tan /3, the blind spots 
occur mostly in the bino/Higgsino region of the plane, in which IM2I < As tan /3 is raised, 
however the bind spots move to lower values of the gaugino mass relative to /i, weakening the 
constraints on the bino/wino parameter space even further. 

Currently, there are no limits from SI direct detection for fi < 0. XENONIT will constrain 
bino/wino DM to lie near the —Mi = /i sin 2/3 blind spot for /i, M2 < 0; for M2 > 0, however, 
the proximity of the well-tempered line to the Mi = M2 blind spot will once again weaken the 
constraints. SD direct detection sets complementary limits, irrespective of the SI blind spots, 
but limits will remain relatively weak in the bino/wino region of the parameter space even after 
XENONIT because of the relatively small mixing angle. 

Regardless of relative signs, the direct detection limits fall off as the magnitude of /i is raised, 
since both the Higgs and Z couplings to dark matter are depleted as /i is decoupled. Even at 
positive n, current limits from XENONIOO only exclude /i < 250 GeV, leaving large allowed 
regions with natural values of without requiring any tuning of the cross-section. LUX and 



30 



current limits 



thermal blw/h 

future reach 




Figure 14: Current limits and future reach for thermal bino/wino dark matter in the {M2, fi) 
plane for tan/3 = 2. Mi is fixed by requiring fl^^ = Qohs- The black dotted lines correspond 
to the SI blind spot for bino-like DM given in Eq. (19). The XENONIOO (XENONIT) SI and 
SD limits (reach) are shown shaded in red and blue, respectively on the left (right) side of the 
figure. IceCube limits on dark matter annihilation to W~^W~ are shown shaded in orange on 
the left, and the LUX SI exclusion reach is shown with a dashed green line in the right-hand 
panel. The LEP chargino exclusion is shaded in light gray. Darker gray regions correspond to 
overclosure via Higgsino-like DM heavier than 1 TeV. 

XENONIT will require /i > 600 and 1500 GeV, respectively, when /i > 0. Evading these limits 
with large fi does not incur a fine-tuning penalty in the cross-section; however, as shown in 
Fig. (15), it is not without cost. Bino/wino DM that evades LUX with fi > must rely heavily 
on coannihilation to achieve the observed relic density. More than 90% (99%) of the total dark 
matter annihilation cross-section, weighted by Boltzmann factors as in Eq. (6), must come from 
coannihilation in order to escape LUX (XENONIT) limits with fi > 0. As discussed above, this 
coannihilation is exponentially sensitive to the mass difference between the bino and the wino, 
leading to a ~ 2% tuning in the relic abundance. The left-hand panel of Fig. (16) shows how 
the relic abundance tuning increases as /i is decoupled. The tuning increases as coannihilation 
becomes more important before reaching a plateau around 1/50. 

By way of contrast, the cross-section tuning, shown in the right-hand panel of Fig. (16), takes 
relatively small values everywhere except close to the blind spots, where it quickly increases to 
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Figure 15: Contours of the fraction of dark matter annihilation cross-section, weighted by 
Boltzmann factors as Eq. (6), coming from coannihilation are shown in purple, superimposed on 
the XENONIT SI exclusion reach shaded in red and the LUX SI reach in dashed green. Once 
again, tan/3 is fixed to 2, with Mi fixed by requiring fl^^ = Qohs- 

~ 1%. Note that in the lower right quadrant, it is possible to evade XENONIT with a relatively 
mild tuning of the cross-section, about 5%, and any value of fi consistent with the LEP bound, 
improving electroweak naturalness. 

7 Conclusions 

In this paper we have presented a systematic analysis of the current limits and projected reach for 
X ~ {b, w, h) neutralino DM using simplified models. For bino/Higgsino and bino/wino (/Higgsino) 
DM we decoupled the heavier supersymmetric partners to yield a parameter space of (Mi, /i, tan (3) 
and (Ml, M2, /i, tan/3), respectively. We assumed CP conservation, but studied all physically 
distinct choices of signs of these parameters. Furthermore, we assumed that the recently dis- 
covered state near 125 GeV is the Higgs boson, and that the other Higgs bosons are decoupled 
and yield subdominant contributions to the DM scattering cross-section. We were then able 
to explore the current limits and future reach directly in this minimal parameter space with- 
out resorting to scatter plots. This simplified model for analyzing neutralino DM is a good 
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Figure 16: Contours of fine-tuning of the relic abundance (SI scattering cross-section) are 
shown in the left (right) panel, using the measure defined in App. B. The XENONIT and LUX 
SI exclusion reaches are shown with solid red and green contours on the right. 

approximation to a wide variety of ultraviolet theories. Limits and reaches were presented for 
neutralinos comprising all or just a fraction of dark matter. The case of thermal freeze-out of 
neutralinos yielding all dark matter was emphasized, as this allows a reduction in the parameter 
space. In what follows, we summarize out main results, and then discuss future directions. 

Our results are detailed in Fig. (3) - Fig. (16), which depict present limits and future reach 
within the theory parameter space of neutralino DM. However, direct detection experiments 
place bounds on the physical {m^, a) plane. For a proper comparison. Fig. (17) depicts the image 
of thermal bino/Higgsino DM in the plane of physical parameters relevant to SI scattering. The 
fi > region has been excluded by XENONIOO for values of up to about 500 GeV, but 
the /i < region is almost entirely unconstrained. The LUX, SuperCDMS and XENONIT 
experiments will probe this fi < region deeply. The absence of a signal would require a 
cancellation in the scattering amplitude at the level of 1 part in 10 - 30. 

For non-thermal bino/Higgsino DM only a small fraction of parameter space with Mi < 
1 TeV has been excluded, as illustrated in Fig. (4). A large (small) fraction of the parameter 
space for thermal DM has been excluded for /x > (/z < 0), as can be seen most clearly in 
the upper panel of Fig. (7). Future experiments, such as LUX and XENONIT, have a large 
discovery potential, as they will explore the majority of the parameter space with /x. Mi up to 
1-2 TeV, as shown in Fig. (5). The hardest region to explore has low tan /3 and /i < 0, as this 
lies close to a blind spot, where Chxx = 0- This is illustrated for thermal DM in the lower panels 
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Figure 17: The gray shaded areas depict target regions in the {m^,asi) plane for thermal 
bino/Higgsino DM, superimposed on the current limit from XENONIOO and the projected 
reaches for LUX and XENONIT. The edge of these gray regions at low results from the LEP 
requirement of \fi\ > 100 GeV, while the largest value of m^, just above 1 TeV, corresponds to 
pure Higgsino LSP, and is present for both signs of fi. The upper dark shaded region is for /i > 
(here we fix Mi > 0) with the upper (lower) edge corresponding to low (high) tan /3. Much of the 
low mass part of this region has been excluded by XENONIOO. The lower two regions, shaded 
in lighter gray, are for /i < 0. The boundary between the fi > and /x < regions occurs at 
large tan /3, where the sign of /i becomes unphysical. In the fi < regions the cross-section falls 
as tan/3 is reduced towards its value at the blind spot, where Mi + sin 2/3 fi = 0. The contour 
between the two /i < regions is given by \Mi + sin 2/3 1 = 0.1 Mi, roughly corresponding to 
a 10% fine-tuning in the scattering amplitude. In the lower region, for each order of magnitude 
further reduction in the cross-section, a factor of \/lO more fine-tuning is required. 

of Fig. (7). Pure Higgsino thermal dark matter will also evade discovery for Mi > 2 TeV, as 
shown by the vertical brown bands in Fig. (5). 

Fig. (6) depicts current limits and projected reaches for bino/Higgsino LSP which is just 
one component of mult i- component DM. Present constraints are quite weak, but LUX and 
XENONIT will probe the fraction of LSP dark matter powerfully, especially at low LSP mass, 
although with the usual blind spot caveat at low tan /3. 

The more general case of bino/wino/Higgsino DM is shown schematically in Fig. (10), and 
contains the interesting possibility of bino/wino thermal DM. Fig. (11) shows the present limits 
and future reach for non-thermal production in a slice of parameter space. While three of the 
four quadrants are affected by blind spots and are currently unconstrained by direct detection, 
all four quadrants will be significantly probed by XENONIT and LUX. Fig. (12) shows the same 
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parameter slice for subdominant thermal DM where, depending on the quadrant, XENONIT 
can probe neutralinos that comprise 10% or even 1% of DM. Fig. (14) shows a slice at tan /3 = 2 
for the present limit and future reach for thermal DM. For /i > (yU < 0) the current limits are 
already quite strong (weak), but even the combination of both SI and SD XENONIT results will 
leave significant regions at yU < with < 1 - 2 TeV. In particular, the large open region in 
the lower right quadrant has Mi close to M2 from the requirement of coannihilation and hence 
lies close to the Mi = M2 blind spot, where Mi 2 have the opposite sign of fi. 

Finally, our results suggest a variety of interesting directions for future work, which we now 
summarize. In particular, 

• The present analysis incorporates the critical effects of relative signs, but still assumes CP 
conservation in the neutralino sector. The importance of these signs for alleviating direct 
detection limits suggest that similar effects can be expected of neutralino DM with arbi- 
trary CP phases. On the other hand, CP-violation in the neutralino sector is stringently 
constrained by electron electric dipole moment (EDM) experiments, subject to the masses 
of the scalar superpartners. It would be interesting to study the interplay between DM 
direct detection and EDM experiments in CP-violating models of neutralino DM. 

• While the bulk of our analysis has been at tree level, our estimates in App. C suggest that 
they are robust to one loop corrections. In particular, naive estimates indicate that the 
projected reach of XENONIT will cover parameter regions for which one loop corrections 
can be likely neglected. However, an explicit calculation of one loop effects in the mixed 
bino/wino/Higgsino system will settle this issue, and have important implications for DM 
in the post-XENONlT era. 

• In this work, we have primarily studied constraints from direct detection experiments and 
neutrino telescopes. An interesting avenue for future work is a more comprehensive analysis 
of the indirect detection constraints on neutralino DM from cosmic ray and gamma ray 
observatories like the Fermi Telescope [29]. 

• Our analysis applies to neutralino DM that is a general mixture of bino, wino, and Higgsino 
DM. However, many of our qualitative results — e.g the importance of relative signs — 
apply much more generally. It would be very interesting to study present and future 
experimental results within the context of generic models of mixed singlet/triplet/doublet 
DM. While simplified models of this type have been considered in a variety of contexts 
[47, 68, 69, 70, 71], their interplay with present and future direct detection limits have been 
less systematically studied; in particular, the analogous blind spot parameter regions have 
not been fully identified. A general singlet /triplet /doublet simplified model would provide 
a theoretically inclusive framework for studying — and, in the absence of positive signals, 
excluding — a huge class of WIMP DM models. Such an analysis would have important 
implications for mixed singlino DM relevant to the NMSSM and its ASUSY variants. 
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It is well known that direct detection limits are sensitive to nuclear physics uncertainties, in par- 
ticular the strange squark content of the nucleon. Higgs-mediated SI scattering is proportional 
to l/p, where [51] 



The largest uncertainty comes from the strange quark content, fs, since fu and fd are small, but 
different determinations have led to widely disparate values for fg, as we now review. Through- 
out, our analysis takes fu = fd = 0.025 [20]. 

Traditionally [72], fg was determined using chiral perturbation theory to relate fs to the 
pion-nucleon sigma term, Ht^n, which is extracted experimentally from the cross-section for 
low-energy pion nucleon scattering, leading to fs = 0.38 ± 0.10. The MicrOMEGAs [51] default 
value, fs = 0.26, is chosen to be near the la lower bound based on these measurements of S^^r. 
However, as first pointed out by Giedt et al. [20], direct lattice determinations of fs lead to 
a significantly smaller value, fs = 0.0532 ± 0.0085. More recent lattice results [73, 74, 75, 76] 
confirm small values, finding fs between 0.009 — 0.046. This ~ 3a tension between the S^tv and 
lattice determinations of fs have led to widely divergent approaches in the theory community, 
with some authors using large values of fs based on S^^v, and others adopting the lattice values. 
We note that this tension is now probably resolved, in favor of the low values, because the most 
recent p — vr scattering data, from the CHAOS group at TRIUMF, lead to fs = 0.10 ± 0.15 [77], 
in agreement with the lattice results. 

In Fig. (18), we show the effect on of the various choices for fs discussed above. The 
cross-section does not change by a large amount between the different lattice determinations. 
However, there is a large difference between the cross-section favored by the lattice and the 
cross-sections resulting from the old determinations of St^at; for example the cross-section is 
increased by a factor of 2.3 when moving between the Giedt and default MicrOMEGAs values. We 
view the lattice determination of fs to be most accurate (especially in light of the new, lower, 
value for S^jv) and throughout this paper we have used the value of Giedt et al. We note that 
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Figure 18: The effect of fs on the Higgs-mediated SI cross-section, which is proportional to 
We show the lattice determinations of fs from Giedt et al. [20] and more recent groups [73, 
74, 75, 76], as well as the value of fs resulting from the most recent measurement of T^^^jy [77]. 
The bands indicate the stated ±1ct errors on the Giedt and E^tat values, which are in agreement 
up to these errors. For comparison, we show the default micrOMEGAs [51] and traditional chiral 
perturbation theory [72] values, which are based on older measurements of Sttat. 

many previous theory studies (for example [19, 78, 17]) have adopted the default MicrOMEGAs 
value, leading to liberal limits by this factor of ~ 2. 

We show the impact of varying fs on the XENONIOO limit, for thermal bino/Higgsino DM, 
in Fig. (19). We compare the Giedt et al. value to the default MicrOMEGAs value and to the most 
conservative possible choice of fs = 0. We see that the choice of fs makes a large difference at 
negative fi and large tan /3, because in this regime the cross-section is close to the XENONIOO 
limit throughout much of the plane. We comment that the thermal bino/Higgsino plane is 
the most sensitive of our results to the precise value of fs, and the limit and reach contours 
throughout the rest of the paper are less sensitive to this choice. 

B Dark Matter Fine- Tuning 

In this paper we are agnostic about the possibility that the EW scale may be finely tuned. But 
whether the weak scale is natural or not, it is interesting to identify when fine-tuning enters 
DM properties. Such a tuning is a worry when DM is well-tempered to produce the correct 
abundance, as discussed in Sec. 3, or if DM sits particularly close to one of the blind spots 
we identified in Sec. 4.2. In this appendix we describe a quantitative measure that tests the 
fine-tuning of fl or a si, independently of a possible EW tuning. We applied the methodology 
described in this appendix to produce figures 8 and 16, above. 



37 




Figure 19: The impact of different choices of fs on the XENONIOO hmit for thermal 
bino/Higgsino DM, which was shown using the Giedt et al. value in the first panel of Fig. (7). 
We compare the Giedt value to the default MicrOMEGAS value and fs = 0, which is the most 
conservative choice. 

To start, let us denote the log quantities, pi which label parameters at the weak scale: 

exppi = {Mi,M2,/i,m^^,m^^,5/i}. (22) 

Then it is natural to define a log gradient defined as a directional directive with respect to log 
parameters, 

V = V, = ^. (23) 

OPi 

We can now define a vector in this space equal to the gradient of the electroweak symmetry 
breaking vacuum expectation value, 

V = Vlogw^ (24) 

Here V is equal to the direction of steepest descent away from a particular value of v — thus, it 
is the combination of ultraviolet parameters which most strongly affects electroweak symmetry 
breaking. In order to remove dependence on the possible fine-tuning of the electroweak symmetry 
breaking sector, we are interested in dependencies on parameters orthogonal to V. Explicitly, 
any v dependent observable can be written as 

= 7^r^^ + ---' (25) 
o log 
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where we have used the chain rule and the eUipses denote dependencies on other parameters. 
Thus our fine-tuning parameter, which is independent of electroweak symmetry breaking, for a 
given observable O is defined by 

Ao = \V^O\, (26) 

where V"*" is defined as the gradient of O within the subspace of the above parameters orthogonal 
to V. In figures 8 and 16, we used A^, where O = ^,<7sh to quantify tuning in the thermal 
relic abundance, and SI cross-section, respectively. 

We compute Ao using tree-level relationships among the weak scale parameters, with one 
exception. The D-flat direction of the tree-level scalar potential causes the aforementioned 
tuning to blow up as tan (3 — )■ 1 (see, for example, [78]). We regard this apparently infinite 
tuning as unphysical, because the D-fiat direction is lifted by loop corrections (and, moreover, 
whatever physics beyond the tree-level MSSM accommodates the Higgs mass). In order to 
lift the D-fiat direction, we assume an MSSM-like completion and add a quartic contribution 
cx: \Hu\'^, with a size fixed to reproduce ruh = 125 GeV. The size of Ao near tan/3 = 1 depends 
on the exact form of potential; however, any completion that lifts the D-fiat direction serves to 
regularize the tuning as well. We have verified that adding a quartic contribution of the form 
iHuHdl"^ yields similar results. 

C Expected Size of Loop Corrections 

To obtain our results we have only included the tree-level interactions of neutralino dark matter. 
In most regions of parameter space this is a good approximation; however, close to the blind 
spots that we identified in Sec. 4.2, the relevant tree-level couplings to the Higgs and Z bosons 
vanish, in which case it is crucial to ascertain the importance of one-loop corrections. While a 
full computation of one-loop corrections is beyond the scope of the present work, the topic of 
this appendix is to estimate when loop corrections would modify our results for SI scattering. 

First, consider the case that DM is a pure eigenstate. If DM is pure bino, it is inert (with 
decoupled scalars) and there are no radiative corrections to the Higgs coupling to the neutralino. 
On the other hand, if DM is a pure Higgsino or pure wino, it has no Higgs coupling at tree-level, 
but there are loop diagrams such as box diagrams with two gauge bosons. The naive size of 
these loop corrections is 10~^^^~^^^ cm^; large enough to probe at XENONIT. However, a dia- 
grammatic calculation including all one-loop and leading two-loop diagrams shows a surprising 
accidental cancellation among scalar and tensor scattering operators generated by one-loop and 
two-loop effects [43]. For example, taking a pure Higgsino (wino) with a mass of 500 GeV, and 
using the form factors of [20], the total cross-section is 6 x 10^^'^ (8 x 10~^^) cm^; much too small 
to probe at XENONIT. Recently, this accidental cancellation was confirmed in an effective field 
theory calculation [44]. 

The multiloop result of Refs. [43, 44] have not yet been generalized to mixed states, as 
would be relevant for the blind spots. We can still estimate the maximum size of the loop 
corrections. Consider the case of mixed bino/Higgsino. The bino component of the DM is given 
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Figure 20: Estimated importance of loop corrections for bino/Higgsino DM, relative to the 
XENONIT reach. Non-thermal and thermal DM are shown in the upper and lower panels, 
respectively, with the XENONIT reach as in Fig. (5) (for tan/3 = 2) and Fig. (7). The green 
contours show the bino fraction of DM, sin^ = Z\y^ and loop corrections are known to be 
small when DM is a pure bino or Higgsino, sin^ Q-^ ?a 0,1. Blue contours show the ratio of our 
estimated loop cross-section to the sum of the tree-level and XENONIT reach cross-sections. 
A large value of this ratio would indicate where the XENONIT reach estimate is sensitive to 
loop corrections, however we see that this ratio is small throughout most of parameter space, 
indicating that our results are robust to loop corrections. 
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by sin^^^ = Zf^, where Zij is the rotation matrix going from interaction to mass eigenstate. If 
sin^ 9i ^ 0, 1, then the DM is close to a pure state and the cross-section is small, as discussed 
above. This quantity is shown by green curves for non-thermal and thermal cosmologies in the 
upper and lower panels of Fig. (20), respectively, and we see that much of the parameter space is 
characterized by a nearly pure state. As the mixing angle is increased, any possible enhancement 
to the cross-section is suppressed at least by a factor of mixing angle squared. Therefore, we 
conservatively estimate the maximum size of the cross-section to be, 

^loop = \Zu\\l - l^iiH X (2 X 10-^^ cm^). (27) 

Here, the number in parentheses reflects the size of the cross-section from the largest individual 
loop diagram contributing to the Higgsino cross-section from Ref. [43] (which happens to be a 
box of W bosons contributing to the tensor operator). We stress that this is simply an estimate 
of the maximum size; the full calculation is beyond our scope. 

When are our results sensitive to the loop corrections? Consider the XENONIT reach, 
shown for bino/Higgsino in Fig. (20). For the loop correction to be relevant, two conditions 
must be satisfied, (1) the loop contribution must be large relative to the tree-level scattering, 
and (2) the loop contribution must be large enough to probe at XENONIT. In order to estimate 
when both of these conditions are met, the blue contours in Fig. (20) show the ratio of crioop 
to the maximum of the tree-level cross-section and the XENONIT limit. The loop correction 
is important when this quantity is large. However, this quantity is less than 0.01 in the entire 
parameter space, except for small regions near the blind spot. Even here, it is only at low fi that 
this ratio reaches 0.05, and its maximum value is near 0.2. Hence, unless our estimate Eq. (27) 
is too small by over an order of magnitude, the loop corrections can be ignored for determining 
the reach of IT detectors for bino/Higgsino DM. 

We have not included an estimate of the importance of loop corrections for mixed DM with 
a large wino component. In this case there are competing effects; the bino/wino DM typically 
has a smaller mixing angle than the bino/Higgsino case, further suppressing any enhancement 
to the loop contribution coming from mixing, but the largest individual loop diagram for pure 
wino scattering is an order of magnitude larger than for the pure Higgsino case. 
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